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ABSTRACT 

Long Gamma-Ray Bursts (GRBs) are produced by ultra-relativistic jets launched from 
core collapse of massive stars. Most massive stars form in binaries and/or in star clus¬ 
ters, which means that there may be a significant external photon field (EPF) around 
the GRB progenitor. We calculate the inverse-Compton scattering of EPF by the hot 
electrons in the GRB jet. Three possible cases of EPF are considered: the progenitor 
is (I) in a massive binary system, (II) surrounded by a Wolf-Rayet-star wind, and 
(III) in a dense star cluster. Typical luminosities of 10 46 — 10 50 erg/s in the 1 - 100 
GeV band are expected, depending on the stellar luminosity, binary separation (I), 
wind mass loss rate (II), stellar number density (III), etc. We calculate the lightcurve 
and spectrum in each case, taking fully into account the equal-arrival time surfaces 
and possible pair-production absorption with the prompt y-rays. Observations can 
put constraints on the existence of such EPFs (and hence on the nature of GRB 
progenitors) and on the radius where the jet internal dissipation process accelerates 
electrons. 

Key words: radiation mechanisms: non-thermal - methods: analytical - gamma-ray 
burst: general. 


1 INTRODUCTION 


Long Gamma-Ray Bursts (GRBs) are produced by ultra- 
relativistic jets launched aft e r th e collapse of mas- 
sive stars fe. g. iGa lama et al.l I 1 9 981: iHiorth et al.l 1 200 .'ll : 

stars 
clus- 


IPiranl 
tend to 


[200 


form 


Wooslev fc BloomT 200a b Massive 

in binaries and/or in dense star 


ters JM assev fc Hunted Il998l: IZ innecke r fc York)] 120071 : 


[Mason et al.l 20091 : Portegies Zwart et al.l l2010ll 7 which pro¬ 
duce a dense external photon field (EPF) around the GRB 
progenitor. Soft photons will be inverse-Compton (IC) scat¬ 
tered by electrons in the jet. The scattering will boost the 
energy of photons by a factor of r 2 7 2 , where P ~ 300 is the 
jet bulk Lorentz factor (LF) and y e is electrons’ LF in the 
comoving frame. Therefore, we expect EPF at 10 eV being 
scattered to 10®7e eV. 

I 11 the GRB literature, this external inverse-Compton 
(EIC) emission has been considered for a variety of sources 
of soft photons and scattering electrons. For example, the 
soft photons may come from a pulsar wind bubble created 
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before the GRB (jGuetta fc Granotli2003i f. and from the fun- 
nel of the collapsing star created during the jet breakout 

jLazzati et, al.l l2000: Ghis ellini et, aLlfeoonrlMacFadven et al.1 


l200lll . Ramirez-Ruid d2004lf considered soft photons from a 
companion star being IC scattered by a magnetically driven 
relativisti c wind from a spinning-down millisecond s upra- 
magnetar. iGiannio s d200£j f ; iMimica fc Gianniosi (1201 il l con¬ 
sidered the GRB external shock interacting with the photon 
field from an O star in the surrounding dense cluster. 


Different from previous studies, we consider three cases 
where significant EPFs exist: the GRB progenitor (I) is in 
a massive binary system, (II) has a strong wind, and (III) 
is in a dense star cluster. 


I 11 Case (I), soft photons come from a massive compan¬ 
ion star. Spectroscopic and direct imaging studies of Galac¬ 
tic massive stars, although biased by selection effects and 
measurement limitations, conclude that massive stars have 
binary fraction of at least ~ 50% (see iMason et al.l |2009| . 
and reference therein). Therefore, a companion of compa¬ 
rable mass may exist near the GRB progenitor, dep ending 
on the evolution status of the binary system flLange ri2012h. 
For binary separation d = 10 15 cm and the companion star’s 
bolometric luminosity Lb = 10 39 erg/s, the number of scat¬ 
tered photons will be 10 52 — 10 53 (isotropic equivalent). 
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In Case (II), soft photons originally from the progen¬ 
itor star are scattered first by the electrons in the stellar 
wind and again by the jet. A compact Wolf-Rayet (W-R) 
star progenitor is favored by considering the p ropa g atio n 
time of a GRB jet in the stellar envelope (M atzner l 120031 1. 
W-R stars have high luminosities of ~ 10 39 erg/s and are 
surrounded by strong stellar win d of a typical mass loss 
rate 10 -5 Mq/uv dCrowtheil2007h . The number of scattered 
photons is 10 5 ° — 10 51 (isotropic equivalent). 

In Case (III), all the stars in the cluster create a 
nearly isotropic and uniform EPF around the_GRB jet. 
Mass ive stars form in clustered environments (lLada fc Ladal 
l2003ll . In our Galaxy, ~ 85% of O stars are either ob¬ 
served in young c lusters or direc tly identified as runaways 
(see IPorteeies Zwart et al.l hold . and reference therein). 
Within the central ~ 0.2 pc of Arches Cluster, the dens¬ 
est known young massive cluster in the Milky Way, the 
number d ensity of massive (30 < M/Mq < 120) stars is ~ 
5000 pc ~ 3 (Espinoza et al.ll2009l ). These stars have luminosi¬ 
ties close to or higher than 10 39 erg/s. Long GRBs are as¬ 
socia t ed with actively star-form ing galaxies ( Fruchter et al.l 
2006; Wo oslev fe Blooinll2006l), where young massive clus - 
ters are particularly abundant (Porteeies Zwart et al.l[2010h . 
Therefore, we expect a good fraction, depending on the clus¬ 
ters’ evolution, of long GRBs to happen in such star clus¬ 
ters. In this Case, we consider EIC emission from both in¬ 
ternal dissipation and external shocks (both reverse shock 
and forw a rd sh ock). Note that our model is different from 
iGiannios! ( 20081 ). who consider the case when the external 
for ward shock happens to swee p across one star in the clus¬ 
ter. iMimica fc Gianniod J201ll) also calculate the EIC scat¬ 
tering of the isotropic diffuse EPF in the cluster, but only 
the external forward shock is considered. As shown in Sec¬ 
tion 3.1.3, the EIC luminosity is dominated by the internal 
dissipation and reverse shock. 

In this paper, we show that all three cases may have 
observable consequences, which can be used to probe mas¬ 
sive stars around GRB progenitors and the nature of GRB 
progenitors. 

A few more notes. (1) By “internal dissipation” (ID) we 
mean the uncert ain process (e.g. interna l shock or magnetic 
reconnection, see lKumar fe Zhandl2015l . for a recent review) 
that dissipates the free energy of the jet and produces the 
prompt 7 -rays. As long as electrons are accelerated to an 
ultra-relativistic powerlaw, our model is independent of the 
details of the jet dissipation and particle acceleration pro¬ 
cess. (2) The opening angle of the jet is assumed to be larger 
than the causally connected cone, i.e. dj t7nax > 1/r, so we 
use the isotropic equivalent energy, luminosity, photon num¬ 
ber, etc throughout the paper. (3) The jet is assumed to be 
hadron dominated and made of pure protons and electron^]. 
(4) The convention X n = A/10 n in CGS units is adopted. 

This paper is organized as follows. In Section 2, we 


1 Hadrons carry the momentum of the jet but are neglected in 
the IC scattering process (since there is an equal amount of elec¬ 
trons). If the jet is magnetically dominated, we expect the EIC 
emission from the internal dissipation to be much weaker, due 
to a much smaller number of electrons. For a multi-element jet, 
e.g. Hydrogen and Helium, our model can be easily modified by 
putting a (1 + X(H))/2 factor before Ej (the total kinetic energy 
of the jet), where X(H) is the mass fraction of Hydrogen. 



Figure 1. Lab frame (progenitor’s rest frame). A spherically 
capped jet is moving in the z direction with bulk Lorentz fac¬ 
tor F. We consider a small volume element (orange) r 2 dfljdr = 
r 2 sindjd9jd(f>jdr of the jet at position r = (r, 0 } , (f>j). The exter¬ 
nal photon field at position r is denoted by intensity I VQ (r. 0, o), 
where (9, <jf>) is photons’ moving direction. Photons scattered into 
the observer’s cone A9 0 f, s are considered as observed. 


give a general treatment of the EIC emission (spectrum and 
lightcurve) when a jet interacts with an arbitrary EPF. In 
Section 3, we first give simple order-of-magnitude estimates 
of the total EIC energy in the three cases, and then put 
the specific EPFs into the general procedures and calculate 
the precise lightcurves and spectra. In Section 4, we discuss 
some potential issues, e.g. absorption of high energy 7 -rays 
by the cosmic background light or by the local EPF itself. 
Short conclusion is given in Section 5. 


2 MODELING THE EIC EMISSION 

In this section, we provide a general calculation for the 
lightcurve and spectrum expected for an observer on the 
jet axis, when a relativistic jet makes its way through an 
arbitrary external photon held (EPF). 

EIC emission has been studied extensively in the 
lite rature, for instan ce in active galactic nuclei (e.g. 
iDermer fe Schlickeiserl 11993T) wh ere a jet me ets a n EPF. 

We refer the readers_to iBlumenthal fe Gouldl (l970l) ; 

lAharonian fc Atoy an|jll981 ) for a general treatment on IC 
radiation. Also, Fan et al.l (2008j) give a review of IC radia¬ 
tion in GRB afterglows, including EIC and synchrotron-self 
Compton (SSC) radiation. Our treatment is different only in 
that, at an early stage, equal-arrival time surfaces are taken 
into consideration. Experienced readers could skip Section 
2 and go to Section 3 for the results. 

In terms of notations, (1) we use un-primed quantities in 
the lab frame (rest frame of the progenitor star) and primed 
(') ones in the jet comoving frame, and quantities with sub¬ 
script “j”, e.g. ( 9j,cf)j ), are related to different regions in 
the jet. Note that, by “comoving frame” we mean “comov¬ 
ing” with the bulk motion of the jet, rather than the host 
galaxy’s redshift. (2) We denote the unscattered and scat¬ 
tered photons’ frequencies as vq (v q) and v (z/) in the lab 
(comoving) frame. 

As shown in Fig. Q. the jet is moving with a Lorentz 
factor (LF) T in the z direction at radius r. For a volume 
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I 

Figure 2. Comoving frame of the volume element (orange) in 

Fig.JTJ. External photon field is denoted by intensity I', (O', ft). 

u o 

The observer’s cone AQ 7 is in the x'z' plane. Photons scattered 
into AO', are considered as observed. 

obs 

element (orange) at position r = ( r,8j,<t>j ), our goal is to 
consider all the EPF that is scattered into the observer’s 
cone Af2 0 i, 3 , which is a very small solid angle in the z di¬ 
rection considering the cosmological distance^]. The EPF 
specific intensity is denoted as 7„ 0 (r, 8,(j>), where r is the 
position of the volume element and (8, <j >) is the direction in 
which photons are moving. 

In the comoving frame of the volume element, the ob¬ 
server’s cone is 


right-hand rule. Therefore, a one-to-one map (8, rf)) —» (8, <j>) 
is obtained from Eq. ©• Then LT gives the intensity in the 
comoving frame 



> = /,„<».«■(£) 

(5) 

and frequency 



vo = 

r(i + pcosd')v ' 0 

(6) 

and direction 



cosd = 

cosd' + p 2 J.' 

1 + PCOS0" * = * 

(7) 


In the comoving ( x'y'z') frame, if electrons are assumed 
to be moving isotropically with LF 7 e , the (average) differ¬ 
ential IC cross section from (vt , 8', </>') to G',e' obs A'o b s = o) 
ifl llAharonian fe AtovanlUosII ) 

d 2 a ^ 3 err r_. z 2 2 z 2 z 2 . 

9v'dfl' oba ~ 16n'y%v' 0 + 2(1 - z) bg( 1 - z) + b 2 e { 1 - z) 2 

(8) 

where z = hv'/{j e m e c 2 ), be = 2(1 — cosd' ic )'y e hv' 0 / (m e c 2 ), 
hv q hv' ^ ■y e m e c 2 be/(1 + be), d' ic is the angle between 

the direction of incident and scattered photons, i.e. cos8' ic = 
sin8' obs cos(/>'sin8' + cos8' ohs cosd', m e is electron mass, c is 
speed of light, and or is Thomson cross section. 

We consider a powerlaw distribution of electrons 

oc 7 f r p d 7 e if 7e ^ 7m 
= 0 otherwise 


A Gobs = D 2 AQ obs 
where Dj is the Doppler factor 


3 r(l — 0cosdj) 


(1) 

( 2 ) 


and p = \J 1 — 1/P 2 is the velocity of the jet. 

Fig.© shows the geometry in the comoving frame, 
where z’ axis is aligned with the original moving direction 
of the volume element. We orient the x axis in a way that 
the observer’s cone is in the x‘z’ plane and at an angle 8' obg 
with the z' axis. Recall that in the lab frame, dj is the angle 
between observer’s cone (same as z axis) and the moving 
direction of the volume element. Therefore, Lorentz trans¬ 
formation (LT) gives 


cos8 obs 


cosdj — P 
1 — p cosdj 


(3) 


Going from 7„ 0 (8, <j >) to 7'/ (8 1 , q /) needs an axis rotation 
and a LT. The axis rotation from xyz frame to xyz frame, 
in which z axis points towards the moving direction, can be 
expressed by 

(el, el, €3) = Ry(8j)R s (n - <j>j)(ei,e 2 ,ez) (4) 

where {el} and {€i} ( i = 1,2,3) are respectively the basis 
vectors in the xyz and xyz frame. Ry,i(ct) are rotation op¬ 
erators, along y or z axis by an angle a according to the 


2 When calculating the isotropic equivalent luminosity, we use 
Af2 0 (, s = 47r. Throughout this paper, if not specially stated, lu¬ 
minosities, emission energies, frequencies and observer’s time are 
presented in the host galaxy frame. 


Therefore, the averaged differential cross section becomes 
f°° „ f'l 7 e ~ p d7 e 


9V 

dv'dQ’, 


= (P~ 1) 


^ym'GGleC 




1 + 


hv' 

1 

2 


POO 

P+1 


7e P d7e 
3 (7t 


( 10 ) 


2 

ae 


16777 

z 2 

1 — z 


m Vq 


+ 


(i - z y 


°dz 


where ae = 2(1 — cos8' ic ) h 2 v' 0 v'/(m e (?) 2 , z ma x = 

vain(hv'/{^ m m e c 2 ), a\ + 4 ae — ae)/2). Note that ae ~ 
h 2 v' 0 v' / {m e c?) 2 is roughly the criterion for whether EIC 
scattering is in the Klein-Nishina (KN) regime. 

In the low energy band v'/p^vb) < 1, the differential 
cross section oc v'° (i.e. constant), so we expect v' L„/ oc v' 2 . 
In the high energy band v'IGlm v 'o) > lj we have 


Zmax 



+ 4 a e - a e 


/2 


1 — 1 /ae + 0{l/a"e) 

\[ae + O(ae) 


if a s » 1 
if ae < 1 


( 11 ) 


Then we get 

/ d 2 o \ f v'~ (p+1)/2 ifa s <l 

\dv'dQ.' ob J^ }v'~ (p+1) {ln{ae) + C) if a e > 1 

( 12 ) 

as pointed out by iBlumenthal fe Gouldl d 19701 1 and 
lAharonian fc Atovanl (jl98ll h Therefore, in the Thomson 


3 The approximation is good if 7 e 1. 
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regime (ag <C 1), we expect u'L„i oc v ,( - 3 p ^ 2 ; in the KN 
regime [ag 2S> 1), we expect v'L u i oc u’ 1 p . However, as 
we shall see in Section the EIC spectrum is almost 
never i/V- p )/ 2 or The reasons are as follows: (1) When 
u /(r 2 lrn 1 ' o) < lj we are generally in the Thomson regime, 
and hence we get vL v oc v 2 . (2) When v / (V 2 ^^ n vo) > 1, 
we almost immediately get into the KN regime {ag > 1), 
and the ( ln(ag ) + C) factor can’t be ignored until ag Js 100 
(i.e. hu ^ 3 TeV ). As the high energy tail suffers from 
pair-production absorption (see the Appendix), the spectra 
hardly reach TeV band. 

Suppose the volume element (orange in Fig[l]) contains 
dN e electrons, the number of photons that are scattered into 
AQ' obs , in a time duration dt' and frequency range du', is 


dN 1 = dN e ■ du' AQ' obs dt' 



du‘. 


,1’AO’A') 


hut 


d 2 


du'dQ' 



(13) 


where the optical depth t(u q, 9') describes the attenuation 
by the part of the jet that lies along the incident photons’ 
trajectory before they enter the jet volume element where 
the IC scattering is considered. Second scattering is ignored 
here. We assume the jet is launched steadily with (isotropic) 
power Lj and duration T), total (isotropic) kinetic energy 
is Ej, and the thickness of the jet is Ar = cTj. Photons are 
all assumed to enter the jet from the front surface, which is 
at radius r. If the volume element is located Sr from the jet 
front, the optical depth along the incident photons’ trajec¬ 
tory is 


r(uo,9') 


Sr Ej , 1 

A rTm p c 2 tot 0 Virr 2 \cos9’\ 


(14) 


where \cos9'\ is due to the incident angle onto the front 
surface and atot{u'a ) is the total KN cross section averaged 
over electrons’ powerlaw distribution 


, iy f-Z aKN ( l '0,7e)7e Pd 7e 

O-tot(Uo) = - 




(15) 


Actually, not all the dN 1 fEg 1131) photons will be ob¬ 
served. Photons of energy hu > 10 GeV may collide with 
low energy photons (such as the GRB prompt emission) 
and produce pairs. In the Appendix, we calculate the pair- 
production optical depth r 7T (^, r) f Eg I A9I) for high energy 
photons produced at radius r and have frequency v. Basi¬ 
cally, we show that only the absorption by the GRB prompt 
emission may be important, if the prompt emission is pro¬ 
duced at t he same radius. And we use a broken powerlaw 
spectrum dBand et al.ll 199^1 with low (high) frequency in¬ 
dex — 1 (—2.4), break frequency 200 keV and bolometric 
isotropic l uminosity 10 5 2 erq/ s , following the Fermi GRB 
statistics (iGruber et al.l 120141) . Then, the diV 7 in Eg. (1131) 
needs to be multiplied by e~ T ’ yl . 

Next we include the effect of the conical shape of the jet 
(“curvature effect”) on the observed luminosity. When the 
jet front is at radius r, the scattered photons from a volume 
element at a distance Sr from the jet front and at latitude 
9j will arrive at observer’s time t obs {r, Sr, 9j). Photons arriv¬ 
ing at the same time are from the same equal-arrival time 
surface, described by 


r eq (6 : j , Sr, fobs) 


0 


ctobs — Sr 
1 — j3cos9j 


(16) 


We divide the jet into many thin shells of thickness d(Sr) ac¬ 
cording to the distance to the jet front Sr and calculate the 
lightcurve contributed by each shell. Using Lorentz trans¬ 
formation in time dt' = dr eg /(Tc) ~ Djdtobs and frequency 
v' = u / Dj in Eq. ( 1131) . we obtain the specific luminosity from 
one shell d(5r) at Sr 


dL , 


_ dN-yhu ^— Tyy ( Uir ) 
dudtobs 

Ej d(5r) f dQ 


rm p c 2 Ar 


huD AQ obs e ' rT 
4"7T j 


-T 77 (i',r) 


(17) 


[ dQ' f 

J4 7T JO 


dun— - 


d 2 a 


hu' 0 \ dv'dQ' n 


obs 
2 t r 




sin9j d9j f 0 n d(j>j and we use 


j :Tnax = 4/r (not sensitive) in our numerical results (Sec- 


where f dQj = f 0 

9 

tion 13.211 . Then, the total observed luminosity is obtained 
by adding up all the shells, i.e. integrating Eq. m over 
f 0 Ar d(Sr). 


3 APPLICATIONS 

In this section, we apply the procedures developed in Section 
[2] to three different cases, where the EPFs are given by the 
physical environment. First, in Section [3.II we give simple 
order-of-magnitude estimates of the total EIC energy, not 
considering the pair production. See Table © for the sum¬ 
mary of our estimates on peak frequency, EIC luminosity, 
duration, etc in all cases. Then in Section m we present 
the EIC lightcurves and spectra, with absorption from pair 
production (with prompt y-rays) considered. For simplicity, 
each case is considered separately in this paper, but in real¬ 
ity, different cases can operate simultaneously (e.g. the GRB 
progenitor is in a binary system and has a wind) and the 
total EIC emission will approximately be a superposition of 
them. 


3.1 Order-of-magnitude Estimate 

3.1.1 Case (I) — Binary System 

Consider the GRB progenitor star being in a binary system. 
Soft photons from the companion star are scattered by the 
electrons in the relativistic jet. 

We assume the direction of the jet is perpendicular to 
the orbital plan«J3- As shown in Fig.®, the jet is at a dis¬ 
tance r from the progenitor star (A) and is moving with LF 
r in the z direction, and the companion star (B) is on the x 
axis at a distance d from A. Photons from Star B encounter 
the jet at an angle cosd * = r/y/r 2 + d 2 . 

Suppose the companion star has a bolometric luminos¬ 
ity L b = 10 39 erg/s L b ,39 and peak energy eo = lOeU eo.i- 
Considering the possible recoil of electrons, we estimate the 
energy of scattered photons by 

e ~ min (V 2 Sl e o, T'y e meC 2 ) 

= 9 GeV r 2 . 5 7 e, 2 eo,i Vkn 


4 Even if the jet direction is inclined by up to 30° from the normal 
direction of orbital plane, the total energy (or luminosity) will be 
affected by a factor of order unity (< 2). 
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Figure 3. Geometry of Case (I). The jet is on the z axis at a 
distance r from the progenitor star (A). The companion star (B) 
is on the x axis at a distance d from Star A. Photons from Star 
B encounter the jet at an angle cosO* = r/y/r' 2, + d 2 . 


Figure 4. Geometry of Case (II). The jet is on the z axis at a 
distance r from the progenitor star. Photons originally from Star 
A are scattered first by electrons in the wind and again by the 
jet. The wind scattered photons encounter the jet with an angle 

e. 


where we have used fiducial LFs T = 300T2.5, 7e = 10 2 7 e , 2 , 
and 


t]kn = min 


17 ) 

r2.57e,2£0,l J 


(19) 


describes the degree to which EIC scattering is KN sup¬ 
pressed. Suppose the jet is launched steadily with (isotropic) 
power Lj = 10 53 erg/s and duration Tj = 10s Tjj , so 
the total (isotropic) kinetic energy is Ej = 10 54 erg Ej : 54 
and the thickness of the jet is Ar = cTj = 3 x 10 11 cm Tj,i. 
Considering possible KN suppression, we estimate the EIC 
optical depth of the whole jet as 


Tj{r) ~ 


Ejar 


4nr 2 rm p c- 


-r/KN 


( 20 ) 


The jet becomes transparent to external photons at radius 


rtr — 3.4 x 10 14 cm 



( 21 ) 


Here, to estimate the total number of EIC scattered pho¬ 
tons, we consider external photons that are in the 7r(9 2 imaa , 
solid angle to be swept by the jet (9j 7max being the half 
opening angle of the jet). As the jet sweeps through the vol¬ 
ume ty9 2 max r 2 dr, the probabiliy of an external photon being 
scattered by the jet is min{ 1, r 2 r /r 2 ). This is because: below 
rtr, the jet is optically thick and all the external photons con¬ 
tribute to the EIC flux; above rtr, the effective cross section 
of the jet is the sum of individual electrons and that equals 
to ^9 2 max r 2 r an d is smaller than 7r 9 2 max r 2 . Note that rtr 
in Ea . (1211) is different from the photospheric radius in that 
the former considers all electrons in the jet but the latter 
only considers electrons in the casualty connected thickness 
r/2T 2 (usually cTj). We note the fact that EIC flux from 
outside the casualty connected region (deeper layers in the 
jet) will arrive at the observer at later time, and this is 
fully taken into account when we calculate the spectra and 
lightcurves (through Eq. 1161) . 

The number density of external soft photons is given by 

= 47r(r 2 + d 2 )e 0 c (22) 


The total (isotropic) number of scattered photons is 

noo 

N-y = / (1 — (3cos9*)n 1 (r, d)min(l, Tj(r))Anr 2 dr (23) 

Jo 

from which, we can easily see that most photons are scat¬ 
tered at radius r ~ d. Therefore, we can estimate 

IV 7 ~ 3 x 10 53 Lb ' 39dl5 min ( 1. % j (24) 

£0,1 V “ / 

Multiplying IV 7 by the scattered photons’ energy from 
Eg. (1181) . we get the total EIC energy 

Eeic ~4x 10 51 erg rl. 57 2 2 Vkn Lb^dismin ^ 1 , ^ 

(25) 

Note that the estimate is not accurate when r tr d, be¬ 
cause the soft photons are moving nearly parallel to jet, and 
hence both N~ t and e decrease by orders of magnitude. 


3.1.2 Case (II) — Stellar Wind 

As suggested by e.g. IWooslevI d 19931 1 and 
iMacFadven fc Wooslevl d 19991 1. GRB progenitor stars 
are possibly Wolf-Rayet (W-R) stars, which have a typi¬ 
cal mass loss rate of 10 M@/yr, bolometric l uminosity 
10 3 9 erg/s, and effective temperature 10 5 K dCrowtheil 
l2007il . We consider that photons originally from the pro¬ 
genitor star are scattered first by the wind and then by the 
jet. The geometry is shown in Fig. ©. 

We assume a steady wind and denote the wind veloc¬ 
ity as v = 10 s cm/s vs and the mass loss rate as M = 
10~ 5 M@/yr M- 5 . The wind-scattering optical depth of the 
region with radius > r is 


Twindir) — 


r 

J r 


M 


4 Tvr 2 vm p c 2 


dr 


~ 2 x iQ-aMziy.- 1 
vs 


( 26 ) 


5 The order of escape velocity. 
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Therefore, the wind-scattered photon density is 

< 27 > 

The number of scattered photons has a logarithmic depen¬ 
dence on r when r < rt r , and then drops off as r -2 at larger 
radius, so we can estimate the total (isotropic) number of 
scattered photons as 

iV 7 = 2.9 x 10 51 Lb • 39? ‘ 1 - 5 In ( (28) 
CO, 1^8 \R*,11 J 

where R* = 10 11 cm R*,u is the radius of the star. There¬ 
fore, the total EIC energy is 


(29) 


f'tr, 14 \ 

R*,u ) 


Eeic = 4.2 x 10 49 er<? 


Lb,39M-5T 2 .5'l e 2VKN , 
- : - In 

vs 


3.1.3 Case (III) — Star Cluster 

Consider the GRB progenitor star being in a young massive 
cluster. All the stars in the cluster create a nearly isotropic 
and uniform EPF around the GRB jet. Scattering by elec¬ 
trons accelerated by both internal dissipation (ID) and ex¬ 
ternal shocks (ES) may be important. We assume that the 
cluster has a(n) (O-)star density of n* = 10 4 pc -3 n t ,4 and 
a radius of R = 1 pc R pc . We assume stars have an average 
bolometric luminosity of Lb = 10 39 erg/s Lb ,39 and peak en¬ 
ergy eo = lOeG eo,i- The number density of EPF can be 
estimated as 

n t L b R 6 -3 n*,4Lb,3gRpc /or .x 

n 7 ~-~ 2 x 10 cm — 1 - : - — (30) 

eoc eo,i 

First, we consider ID. After the jet becoming transpar¬ 
ent (r > rtr), 4-7T r 2 n~/Tj is a constant. The total number of 
scattered photons depends on the largest radius where elec- 
trons stay hot r max , which is quite u ncertain (IPira nl l2004 
iMeszaro a l2006l : iKumar fc Zhandl2015l b Observationally, the 
duration of GRB prompt emission (Tqq/(1 + z) ~ 10 s) gives 
an upper limit of the radius where the ID is active 

rio < 2r 2 cT 9 o/(l + z)~ 5 x 10 16 cm r 2 2 . 5 (31) 

Theoretically, an upper limit of rm is the deceleration ra¬ 
dius rdec (see Ea l37l below'l. where the jet starts to decelerate 
substantially. Note that, after the prompt emission is pro¬ 
duced, electrons cool via adiabatic expansion 7 e oc r -2//3 
and may stay hot for some time. Putting the uncertainties 
in the parameter r ma x = 10 16 cm r ma x,i6, we estimate the 
total EIC energy by 

Eeic —{erjU^^nr 3 /3)\r=r ma x 

-1-4 x 10 5 O er<j E jM T 2 .s'll, 2 , nKN (32) 

’hi* ,4-^6,39 Rpc T max, 16 

Next, we consider the EIC emission from ES, where 
both the reverse shocked (RS) and forward shocked (FS) 
region may contribute significantly. EIC emission from ES 
strongly depends on the ES dynamical evolution. 

We assume pressure equilibriuir@ between FS and RS 


6 This may only hold for a so-called short-lived RS, but see lUhml 
1120111 1 for a description of a long-lived RS, where RS may reach a 
much larger radius. Then the EIC energy from RS could be much 
larger. 


region (lKobavashil2000l l . both moving together with LF r s h- 
Then V s h is only a function of jet LF T and the density 
ratio n/n e j, where n is the number density of the circum- 
burst medium and n e j is the comoving number density of 
unshocked ejecta. Before RS crosses the jet, the LF of the 
shocked region is 

r ‘ h(r) = (1 + 2V^)K~K)Y^ (33) 


The comoving number density of the unshocked ejecta is 


n e j (r) 


4,nr 2 TjF 2 m p c 3 


(34) 


The circum-burst medium is either a uniform density 
medium or stratified like a wind (n oc r -2 ) out to the 
wind termination shock radius. GRB afterglow analyses sug¬ 
gest that about half of the long GRBs have a uniform 
or weak (1/100 W-R) wind de nsity profile, ra ther than a 
normal W-R wind p rofile (e.g. IPanaitescu &i Kumail l200ll ; 
iChevalier et alj|2004f ). the reason for which is still an open 
question. 

For simplicity, we take the electrons’ number density in 
the cluster to be uniform n = 10cm -3 m. Then the jet dy- 
na mics is chara cterized by the deceleration radius rdec (see 
e.g. IPiranl [20041. At radius r < rdec, whether RS is relativis¬ 
tic (r s h < r) or Newtonian (r 3 h « T) depends on the factor 
in Eq.® 2P {n/n e j) 1/2 ~ QASrlri 7 {niTj,i/E jM f ' 2 . For a 
fast jet, say T ^ 10 3 , from Eg. (13311 . we have 


r s fe(r) ~ 1.4 x 10 2 r 17 1/2 


f Ej, 54 \ 

Ui ^-,1 / 


1/4 


(«r) 


(35) 


For a slow jet, say V = 100, 2T(n/n e j) 1 ^ 2 ^ 1, so we have 
r s ft(r < r d ec ) ^ 100. 

The deceleration radius can be estimated by 

r 2 sh ^-r 3 dec nm p c 2 = Ej/ 2 (36) 

which gives 


E 1/3 

r dec cs 9.3 x 10 16 cm ^\ /3 

sh,2 n l 


(37) 


Therefore, for a jet LF V £ (10 2 ,10 3 ), we always have 
r d ec ~ 10 17 cm and r dec = r s h(r dec ) ~ 100. Hereafter we 
use r d ec = 10 17 cm r de c, 17 and F d ec = 10 2 Pdec ,2 as the fidu¬ 
cial values. 

Electrons in the RS and FS region are accelerated to 
different LFs. Assuming a fraction e e of the shocked fluid’s 
internal energy goes into electrons, we have 



(38) 


where m p and m e are proton and electron mass. Below, we 
consider RS and FS separately. 

(1) RS: EIC emission from RS peaks when RS crosses 
the end of the jet. After shock crossing, a rarefaction wave 
propagates through the RS region at sound speed and de¬ 
creases the internal energy substantially. Despite the dif¬ 
ference between shock crossing time (slightly earlier, typi¬ 
cally) and deceleration time (later), a good approximation is 
that RS region stays hot until r de c■ The energy equipartition 
parameter e e , r is poorly constrained from GRB afterglows, 
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Table 1. Summary of the EIC emission from our analytical calculation (Section 13.11 not considering pair-production absorption and 
electrons having a single LF), including hv p (peak energy), (isotropic EIC luminosity), t Q b s (duration of EIC luminosity), 

IV eic (i so f ro pi c number of EIC photons at ~ v p ). The scaled parameters are all set to 1, including Ej y 54 (isotropic jet kinetic 
energy), Tj t i (jet duration), 1 ^ 2.5 (jet bulk LF), id 2 (ID-accelerated electrons’ LF), r ma;C) i 6 (the maximum radius electrons stay 
hot), T^ec 2 (bulk LFs of FS and RS at deceleration radius), 7 e ns 2.5 (RS-accelerated electrons’ LF), r^ ec 17 (deceleration radius), T 4.7 
(temperature of EPF-contributing stars), L 5 39 (bolometric luminosity of EPF-contributing stars), M_5 (wind mass loss rate), vs (wind 
velocity), 4 (stellar number density in the cluster), R pc (radius of the cluster). 


Cases 

o 

II 

/ (d - 10 15 ) 

I (dm 10 16 ) 

II 

III (ID) 

III (RS) 

III (FS) 

hu p [GeV] 

9 

9 

9 

9 

9 

9 

26 

l e°ic i er 9/s] 

4 x 10 49 

5 x 10 49 

5 x 10 48 

5 x 10 48 

10 49 

5 x 10 48 

2 x 10 46 

tobs [*s] 

10 

10 

10 

10 

10 

2 x 10 2 

2 x 10 4 

\riso a 

3 x 10 52 

3 x 10 52 

3 x 10 51 

3 x 10 51 

10 52 

10 53 

6 x 10 51 

r ef fl \ cm \ 

d(10 14 ) 

d(10 15 ) 

d(10 16 ) 

r tr (3 x 10 14 ) 

rmai(10 16 ) 

rdec(10 17 ) 

r KN ( 3-3 x 10 17 ) 


a 10 52 photons from 1 Gpc (z = 0.25) give a photon fluence of 10“ 4 cm -2 . With a collecting area ~ 10 4 cm 2 , Fermi Large Area 
Telescope (LAT) is capable of probing the EPF around low-redshift (z < 0.5) GRB progenitor stars. Also note that he EIC photons 
from the forward shock in Case (III) span a large range of energies. The isotropic number of EIC photons in 0.1 — 300 GeV (Fermi LAT 
sensitivity range) is ~ 2.5 x 10 53 and the duration is t 0 b s ~ 9 d (see Section 13.1.31) . 

^ By using a constant 7 e at all radius in our analytical and numerical calculations, we are making an assumption that electrons are hot 
near the radius r e ff where most scatterings happen (see Fig[6]for a justification). 


since there are not many confirmed RS detections. Consider¬ 
ing the uncertainties, we use e e , r = 0.1 (the same as in FS). 
Therefore, from Ea. (f38|) . we get 7 *. ~ 1.8 x 10 2 (1 + T/2T s h), 
i.e. typically a few hundred. Hereafter, we use 7 e = 3007 e , 2.5 
as the fiducial value for RS. 

According to Eq. G9, the scattered photons’ energy is 
e ~ 9GeV r| ec 2 7 e, 2 . 560,1 (marginally KN suppressed) and 
the total (isotropic) EIC energy from RS can be estimated 
as 

Eeic = (eTjn^nr 3 /3)\ r =r dec 

- 1-4 x 10 51 erg ^P^-n tA L bt3a R pc T 2 d 2 ^ 2Ji r dec,i7 ^ ^ 

i 2.5 

(2) FS: We use Eq. (1381) with a constant e e = 0.1 
dPanaitescu fc Kumar! 120011 1 to get the electrons’ LF 7 e . 
At rdec, EIC radiation from FS region is strongly KN sup¬ 
pressed, so the total EIC energy is rising with r (Eeic c* 7-4 
as shown below), until the scattering changes from KN to 
Thomson regime. The characteristic radius where this tran¬ 
sition happens is denoted as tkn, given by 

TshirKN^eto = m e c (40) 

After r^ec, we expect t he Blandford-Mckee evolution 

(jBlandford fc McKeelll976l 'l 

/ \ ~ 3 / 2 

T sh (r) ~ r dec I — (41) 

\ f'dec / 

Putting Eq. (l38l) and (|4TI) to Eq. (l40l) , we get the “KN radius” 

r KN ^ 3.fr dec r‘£, a #* (42) 

We know E^fc x er F sr 3 . When rdec < r < rKN, the 
scattered photons’ energy e oc r~ 3 , and the optical depth of 
FS region tfs oc r 4 , so we get Eeic x 7,4 • In a similar way, 
when r > rKN, we get Eeic 7,-2 • Therefore, most EIC 


energy is produced at radius rKN, and we obtain 

Eeic — (eT F sn-, Yur 3 /3)\ r=r KN 

in50 n*,4,Lb,39Rpc 4 ^8/3 (43) 

- 2.6 x 10 erg ’ " nir dec ,i 7 r/ 2 

e o;i 


where we have used e(rxiv) — (m e c 2 ) 2 /eo = 26GeE , and 
tfs(vkn) — rKNncjT /3. We note the strong dependence on 
both rdec and Tdec■ Also, higher energy (up to ~ TeV) EIC 
photons are produced at radius r dec < r < rKN, but they 
contribute a smaller EIC energy. At radius r > rKN, the 
decreasing of EIC energy Eeic ^ r ~ 2 is mostly due to the 
fast decreasing of EIC photons’ energy e oc T 4 oc r^ 6 . On 
the other hand, the total isotropic number of EIC photons 
is rising with radius Neic = ^r 3 T F sn^/3 oc r 4 and so 
as the observational duration t 0 b s oc r/Y 2 oc r 4 . Therefore, 
if we use the the same fiducial parameters as in Eg. (1431) . 
( 1 ) the total isotropic number of photons in the 0.1 — 
300 GeV energy range (Fermi Large Area Telescope sensi¬ 
tivity) is Vg°c( 7 \KTA r )(26/0.1)‘P 6 ~ 2.5 x 10 s3 ; (2) the obser¬ 
vational duration is tobs — (26/0.1) 4/,6 rKjv/[2F 2 h (rxAr)c] ~ 
9 d and the photon num ber flux is constant with time. 
iMimica fc Giannios! (l201llf also calculated the EIC scatter¬ 
ing of diffuse EPF in a dense star cluster by the forward 
shock. They used a lower stellar number density and a more 
realistic circum-burst medium density profile (by simulating 
the wind-wind collision between a Wolf-Rayet star and an 
O star), and the results in this paper are consistent with 
theirs. 


3.2 Lightcurves and Spectra 

In this subsection, we put realistic EPF profiles into 
the procedure developed in Section [2] and calculate the 
EIC lightcurves and spectra for the three cases. As a 
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Figure 5. Spectra at t Q b s = 3 s (upper panel) and lightcurves 
at Up in the Binary Case, for different binary separations d and 
bulk Lorentz factors F (e.g. “dl4” means d = 10 14 cm and “g2.5” 
means F = 10 25 = 300). (v p s are the peaks of vL v spectra with 
absorption, solid lines, in the upper panel.) We show the dif¬ 
ferences between calculations with/without pair-production ab¬ 
sorption by using dashed/solid lines (“_a” means absorption is 
included). Absorption is weaker for a larger bulk Lorentz fac¬ 
tor or a larger binary separation. The scaled parameters Ej 54 , 
1 , 7 m, 2 , T 4 . 7 , Lb ,39 are all set to 1 and the electron power- 
law index p = 2.2. Th e thic k line is the “classical” Band 

spectrum (|Band et al.ll 1993 ) with E pea & = 500 keV,ot = —1,/3 = 
—2.4, Ltot = 10 52 erg/s. The thick dashed black lines are differ¬ 
ential sensitivity of different telescopes. The EIC spectra show 
uL v oc v 2 below v p and oc i /~~ 0,5 above v p , plus additional high 
energy softening caused by pair production. 


general assumption, in each case, the star(s) that con¬ 
tribute to EPF are assumed to have bolometric luminos¬ 
ity Lb = 10 39 erg/s Lb ,39 and effective temperature T = 
5 x 10 4 K T 4 . 7 , so the specific luminosity is 


u ° csbT 4 


tv Lb 2hvQ 1 


q 2 ghvQ/kT _ 1 


(44) 


where usb is the Stefan-Boltzmann constant. 

(I) Binary Case : We consider the companion star as a 
point source, so the EPF intensity at position f = (r,9j,ij>j), 
in the direction (9,4>), is 


K 0 {r, 9, <j>) = sin9*5{9 - 0»)<5(0 - n) (45) 


where sind * = d/y/r 2 + d 2 , d is the binary separation, and 
5(x) is the Dirac-5 function. 

(II) Wind Case : Assuming a steady wind with mass loss 
rate M and speed v, we get the number density of circum- 
stellar electrons at a distance r from the center 


n e (r) 


M 

4nr 2 v 


(46) 


Figure 6. Spectra and lightcurves for different jet parameters 
in the Binary Case. In the fiducial case (red), Ej 54 , Tj t 1 , T 2 . 5 , 
7 m, 2 , ^ 15 , T 4 . 7 , Lb 39 are all set to 1 and the electron power- 
law index p = 2.2. Magenta lines are for 7 m = 10, so the scat¬ 
tered photons have lower energies and the un-absorbed spectrum 
is slightly rising from 0.3 to ~ 10 GeV, because the scattering 
is in the Thomson regime (uL u oc p )/ 2 ). Green lines are for 
7 m = 10 3 , so the scattered photons reach ~ 300 GeV, but ab¬ 
sorption is strong. Blue lines correspond to p = 2.8, and the only 
difference from the red lines is that the high energy spectral slope 
is steeper, approaching u 1 ~ p faster because the /n(a^) term (Eq. 
112D plays less of a role. For cyai lines, we assume electrons are hot 
only in a limited radius range r E (d/2, 2d). The small difference 
between the red and cyan lines means that, by using a constant 
7 m at all radius, we are making an assumption that electrons are 
hot near the radius where most scatterings happen (see r e ff in 
Table [lj. The line in the upper panel is the “classical” 

Band spectrum (see Fig[5]for details) shown for comparison. 


Then the emissivity from wind scattering is 


>oM 


n e (r)a T 

4nr 2 


■ L 


"0 


(47) 


which is assumed to be isotropic (no dependence on 6 or <j>) 
from Thomson scattering. The system is spherical symmet¬ 
ric, so jv 0 doesn’t depend on 8j or (f>j , and hence neither does 
the intensity We integrate along the ( 9 , <j>) direction and 
get the EPF intensity at position f 

POO 

Ivo( r ,9, (j>) = / j v (r(s))ds 

Jo (48) 

-»n i n 9 ^ - 9 + sinOcosQ L uo 

vs sin 3 9 r 3 

which is independent of tj>, due to symmetry. Note that the 
function f(9) = ( 7 v — 9 + sindcosd)/sin 3 9 ~ 2(7r — 9) 2 / 3 —¥ 


































Probing Massive Stars around Gamma-Ray Burst Progenitors 9 


0 when 6 —> n~ and f(6) ~ 9~ 3 —> oo (diverge^) when 
9 —y 0 + . However, photons with 9 < 1/T are moving nearly 
parallel to the jet and hence don’t contribute significantly 
to the spectra. 

(Ill) Star Cluster Case: All the stars in the cluster cre¬ 
ate a more or less uniform and isotropic EPF, so the intensity 
can be estimated as 


TI*Li/qRc 

- 47T 


(49) 


where n» is the number density of stars and R c is the radius 
of the cluster. 

In Fig. <[5j)-([ 8 ]), we present the lightcurves and spectra for 
the three cases, assuming electrons are accelerated by inter¬ 
nal dissipation. Then we discuss the prospects of detecting 
the predicted high energy photons. We show the represen¬ 
tative spectra (upper panels) at observer’s time t 0 b s = 3 s 
for Case (I), t 0 b s = 1 s for Case (II), and t 0 b s = 10 s for 
Case (III)- Since the overall flux levels of the spectra change 
with time (but the shapes change very little), the different 
tobs’s we have chosen are when the spectra have their average 
amplitudes. By integrating the representative spectra over 
the duration of EIC emission ~ T, = 10 s, we obtain the 
total fluence approximately. The jet parameters Ej t 54 = 1, 
Tj : i = 1 , T 2.5 = 1 , 7 m ,2 = 1 and p = 2.2 are our fiducial 
values in all three cases. To explore the effect of varying jet 
parameters on EIC emission, we also show the results for a 
number of different parameters (one at a time): T = 1000 
(P 3 = 1), 7m = 10 (7m,1 = 1 ), 7m = 1000 (7m,3 = 1) 
and p = 2.8. For each set of parameters, the solid (dashed) 
line is when pair-production absorption with GRB prompt 
emission is (not) considered. Then we measure the peak fre¬ 
quency v p of the solid line (with absorption) for each set of 
parameters. The lightcurves (shown in lower panels) are cal¬ 
culated at u p . In this way, the solid-line lightcurves are close 
to bolometric and we can see how strong the absorption is 
at Up. 

To demonstrate the detectability of the IC emission, 
(1) along with the EIC spectra, we also plot the “classical ” 
prompt emission, i.e. Band spectrum (iBand et alj 1 19931) . 
with Epeak = 500 ke V, a = —1 ,f} = —2.4, L to t = 10 52 erg/s. 
It’s still an open question whether the Band spectrum ex¬ 
tends to a few GeV without a cutoff. In some cases, the EIC 
emission is only observable if Band spectrum has a GeV 
cutoff. (2) We show in Fig. ©-© the sensitivit y of Fermi 
Large Area Telescope (LAT. lAtwood et al .112009 1). High Al - 
titude Water Cherenkov iHAWC. lAbevsekara et al.l 2012j~) 
obser vatory, Cherenkov Telescope Array fCTA. Ilnoue et all 
l201fj ) with integration time 10 s for a GRB at redshift z = 
0.34, whic h is the redshift of the n earest LAT detected GRB 
(130427a, lAckermann et ai1l2014T ). 

The sensitivity of Fermi LAT is calculated using 
uF v \ min = huN ph /(A e ff(u) * Untg), where N ph = 1 is the 
minimum number of photons, ti n t g = 10 s is the integra¬ 
tion time, and A e ff(u) is the detector’s effective area (for 


7 The divergence is caused by the assumption of the star being a 
point source. If we denote the progenitor star’s radius as IU. at 
an angle 9 < arcsin(R+/r), photons come directly from the star. 
However, as long as the star is not a giant (it* < 10 13 cm), the 
photons directly from the star are moving nearly parallel to the 
jet and are hence negligible. 




Figure 7. Spectra and lightcurves for different jet parameters in 
the Wind Case. In the fiducial case (red), Ej } 54 , Tj 1 , IVs, 7 m, 2 , 
T 4 . 7 , f. 7 , 39 , M— 5 , vg are all set to 1 and the electron powerlaw 
index p = 2.2. Magenta lines correspond to 7 m = 10, and the 
scattered photons have lower energies. Green lines are for = 
10 3 , so the scattered photons reach ~ 300 GeV but the absorption 
is strong. Blue lines correspond to p = 2.8, and the only difference 
from the red lines is that the high energy spectral slope is steeper, 
approaching v faster because the ln(ag) term (Eg. 1121) plays 
less of a role. lines are for T = 1000. The line in 

the upper panel is the “classical” Band spectrum (see Fig[5]for 
details) shown for comparison. 


“sour ce” class photons at normal incidence, [Atwood et al.l 
120091) . Considering background and non-normal incidence 
angles, the true sensitivity of Fermi LAT is worse than 
shown here. We can see that Fermi LAT doesn’t have enough 
collecting area to effectively constrain the EIC emission pa¬ 
rameters for GRB with a > 0.34. However, HAWC and CTA 
will be sensitive enough. 

lAbevsekara et akl (120151) calculate the HAWC (-300) 
sensitivity with integration time tintg = 21.5 s, and we scale 
their results for the main data acquisition system by t; n 7 9 
and obtain the sensitivity for 10 s intetration ti me. As f or the 
CTA sensitivity, we follow the procedure of lBernlohr et all 
(l2013h and calculate the minimum flux level of CTA South 
that satisfies the following three conditions: significance level 
a ^ 5, there are at least 10 signal events, and the signal is 
at least 5% of the remaining background in the direction of 
the source. 

A summary of our numerical results is as follows: 

(1) EIC emission is concentrated in 1 - 100 GeV band. 
The EIC energy is typically 10 48 ~ 10 51 erg (hence lu- 
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Figure 8. Spectra and lightcurves for different jet parameters 
in the Star Cluster Case. In the fiducial case (red), Ej 54 , Tj 1 , 
r 2 .5, 7m,2, ^4.7) 39i r ma x, 16; 71*, 4 , Rpc are all set to 1 and 

the electron powerlaw index p = 2.2. Magenta lines correspond 
to 7 m = 10 , so the scattered photons have lower energies and the 
un-absorbed spectrum is slightly rising from 0.3 to ~ 10 GeV, 
because the scattering is in the Thomson regime. Green lines are 
for 7 m = 10 3 , so the scattered photons reach ~ 300 GeV, but 
the absorption is strong. Blue lines correspond to p — 2.8, and 
the only difference from the red lines is that the high energy 
spectral slope is steeper, approaching v 1 " p faster because the 
ln{ag) term (Eg. 1121) plays less of a role. Cyai lines are for T = 
1000. For black lines, the maximum radius that electrons stay 
hot is r max = 10 15 cm. The line in the upper panel is 

the “classical” Band spectrum (see Fig[5] for details) shown for 
comparison. 


minosity 10 4 ' — 10 5 ° erg/s ) and the photon fluence is 
10 -7 ~ 10 -4 cm ~ 2 for a typical redshifl[f| z = 1. 

(2) The numerical spectra follow the scaling of the av¬ 
eraged differential cross section (Eg 1101) quite well. Without 
considering pair production, we get vL v oc v 2 below v p and 
oc v~~ 0 ' 5 above u p . In the cases where 7 m = 10, a short seg¬ 
ment of vL v oc v^-v)/ 2 can be seen, because the scattering 
is in the Thomson regime. At higher energy, the spectra be¬ 
come increasingly steeper and the expected z/ 1_p slope only 
shows at ~ TeV. However, pair production usually makes 
the high energy tail much softer. 

(3) The observational challenge is that the EIC emis¬ 
sion from the internal dissipation arrives within ~ 10 s 
of the GRB trigger, which requires that the GRB occurs 
in the large field of view of a telescope, e.g. Fermi LAT, 


8 The luminosity distance D ^,. = i = 2.0 X 10 28 cm from Dl = 
c(l + z)/H 0 J 0 Z dx/ \/f2m(l + x) 3 + Ha, and U m = 0.27, 12a = 
0.73, H 0 = 71 km. s~ 1 Mpc~ 1 . 


HAWC and wide-field-mode CTA. However, in the Star 
Cluster Case, although not shown in the figures, EIC pho¬ 
tons from the reverse shock will arrive at deceleration time 
tdec = r de c/(2cr 2 dec ) ~ 200 s and those from the forward 
shock will arrive at the time when the scattering changes 
from KN to Thomson regime tKN = rKN/[2cT a h{rKN) 2 ] ~ 
2 x 10 4 s. This gives enou gh time f or some telescopes , 
e.g. CTA dlnoue_eLaL 2013 ), MA GIC (I Albert e t al.ll2007l). 
H.E.S.S. (lAharonian et alT 20091) .VERITAS ( Acciari et al l 
B)> to carry out follow-up observations. 

(4) In the Binary Case, the strongest EIC emission is 
expected when d ~ 10 14 — 10 16 cm.. Therefore, observations 
can probe the existence and properties of a possible lumi¬ 
nous companion in this distance interval, as long as electrons 
are hot at radius r ~ d. M any GRBs have been d etected in 
the 0.1 — 300 GeV band (lAckermann et aTI 1 20131) . but the 
radiation mechanism of the GeV photons durin g the ea rly 
prom pt emission phase is still uncertain (iKumar fc Zhana 
120151) . so the observed data points can be used as upper 
limits unless different components can be separated and 
subtracted. Among all the low re dshift GRBs with LAT 
GeV data, GRB 130427a (z = 0.34. lAckermann et al.ll2014l ) 
provides the most stringent constraints. However, for the 
most favorable jet parameters, the quantum limit of LAT 
is still slightly above the predicted fluxes. Therefore, the 
constraints from LAT are weak. For example, assuming jet 
parameters and binary separations as in red, magenta and 
blue lines in Fig. [5] we can rule out a companion star with 
luminosity > 5 x 10 39 erg s _1 . On the other hand, HAWC 
observatory should have been able to put strong constraints 
on the external photon field around the progenitor of GRB 
130427a. Unfortunately, the burst occurred under unfavor- 
able observational cond itions (much smaller effective area, 
lAbevsekara et aI1l2015l) . For example, non-detections under 
normal observational conditions would rule out a companion 
star of luminosity > 3 x 10 38 erg s _1 at d £ (10 14 , 10 15 cm), 
if T 500 and 7 m ~ 100. 

We note that, in a possible binary Population III (Pop 
III) star system, the lum inosity of the companion star could 
be as high as 10 4 ° erg/s dHeger fc Woosle^l2002l) . A Pop III 
GRB could also have a higher jet kinetic energy Ej : 54 = 10 
(ISuwa fc Iokall20lH ). If a GRB happens in a Pop III binary 
systen@ where the separation between stars is d = 10 15 cm, 
for our fiducial parameters, we expect N^ic — 2 x 10 53 at 
173/(1 + 2) ~ 8 /(l + 2) GeV, which means a photon fluence of 
~ 10 ' cm -2 at redshift 2 = 15 and that could be observed 
by future large space telescopes. 

(5) In the Wind Case, EIC emission is usually weak 
( Eeic ~ 10 49 erg), and observations by LAT are not very 
constraining. Since Eeic oc (Eg 1291) . this case is worth 
considering because of two possible enhancement channels. 
First, there might be a persistent (months to years) and 
strong ( 10 -3 ~ 1 Mg/yr) pre-burst mass loss, as seen in 
some Type Ibc and all Type Iln supernovae (SNe, e.g. 
iFolev et aPl2007l : ISmartti 2009; Kiewe et al.ll 2012 l) . We also 


9 The multiplicity and binary separation distribution of Pop III 
stars are largely unknown. By simulating the growth and evolu- 
tion of Pop III stellar s ystems in a sample (N = 10) of minihalos, 
IStacv &; Bromml (l2013h find a binary fraction of ~ 35% (at birth) 
and that the binary separation (at birth) peaks at 10 2 AU. The 
subsequent evolution of these binary systems is still unexplored. 
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note that the strong mass loss might not be isotropij 19 ! and 
hence could leave no footprints on the afterglow dynamics. 
Second, an unusually high stellar luminosity might be pro¬ 
duced by, e.g. a SN slightly earlier than the GRB, or a hot 
cocoon 11 ! surrounding the jet. 

( 6 ) In the Star Cluster Case, the EIC emission from 
internal dissipation is weaker than the Binary Case by a 
factor of ~ 10. We emphasize the importance of EIC ra¬ 
diation from the external shock (ES), because we have a 
better understanding of the ES than of the internal dissi¬ 
pation in that the standard (low energy) GRB afterglows 
are most likely due to syn chrotron ra diation from the ES- 
accelerated electrons (e.g. IPirani [2004 1. However, the un¬ 
certainties of the EIC emission from ES come from the 
jet dynamical evolution, which depends on circum-burst 
number density and jet power. Also, GeV emission from 
other mechanisms needs to be separated, e.g. synchrotron 
and synchrotron self-Comp ton emission from forward shock 
(see iKumar fc Zhanj|2015l . for a review). We predict pho¬ 
ton fluence from 10 ~ 6 crrU 2 (forward shock, ~ 30 GeV) to 
10 -5 cm~ 2 (reverse shock, ~ 10 GeV) at redshift z = 1, us¬ 
ing our fiducial parameters. Future large telescopes may be 
able to constrain the properties of the cluster. 


4 DISCUSSION 

In this section, we discuss some potential issues of the EIC 
emission proposed in this paper. 

(1) We note that high energy 7 -rays can’t propagate 
through cosmological distances, due to pair production with 
extragalactic background light (consisting mostly of cosmic 
infrared, optical and UV backgrounds). Photons of energy 
(S 100 GeV can onl y be observed from nearby (z < 1) GRBs 
JCilmore et al.ll2012l l. Therefore, constraints on the external 
photon field will be most effective for low redshift (z < 1 ) 
GRBs. We also note that GRBs with low bulk Lorentz factor 
r < 300 will mostly produce < 10 GeV photons, which could 
be observed at high redshift. 

(2) We assume the jet kinetic power to be steady for a 
duration Tj. In reality, GRB prompt lightcurves show mul¬ 
tiple emission pulses. Some pulses are smooth but others 
vary on millisecond timescale, so the jet kinetic power could 
be episodic and the internal dissipation could happen lo¬ 
cally. The EIC emission is also modulated by the radial dis¬ 
tribution of external photon held. If most EIC scatterings 
happen at radii ~ r, an infinitely thin shell will produce 
a pulse of width ~ 0.2s ns/rl.s- Due to limited number 
of EIC photons from a source located at cosmological dis¬ 
tances, telescopes are not sensitive to the variability down 
to milliseconds. Instead, the observed EIC lightcurve should 
more or less track the major pulses of the prompt emission, 
with smoothing on timescales of order ~ 0 . 2 s ris/F 2 5 . 

(3) In the Binary Case, the EIC luminosity is only sig¬ 
nificant when the binary separation is large (d > 10 14 cm). 
Orbital separation distribution derived from spectroscopic 


10 For example, due to the fast rotation of the progenitor, most 
of the gas coul d be ejected nea r the equatorial plane. 

11 See lKumar fe Smoot! ||2014| 1 for the EIC scattering of the hot 
cocoon radiation by the GRB jet. 


and direct imaging studies are inevitably biased by selec¬ 
tion effects an d measurement li mitations. By measuring ra¬ 
dial velocities. ISana et alj l)2012| j derived (after completeness 
corrections) the distribution of orbital periods f(logP) oc 
( log P) -0 ' 55 f or p = 1.4 — 3200 d by analyze the O star pop¬ 
ulation of six nearby G alactic open clusters. Al so by mea¬ 
suring radial velocities, iKobulnickv et all ([2014) estimated 
f{log P) oc ( log P) -0 ' 22 f or p = 1.4 — 2000 d in the Cygnus 
OB2 Association. Despite the discrepancy (and also uncer¬ 
tainties), both of them pointed out that about 10 % of the 
binary systems have periods P > 10 3 d, which corresponds 
to d ~ 10 14 cm. An even greater uncertainty comes from the 
unkno wn nature of G RB progenitors, which may n ot be O 
stars llWooslevI Il993l : iMacFadven fe Wooslevl 1 1999! ) . Detec¬ 
tion of the EIC emission proposed in this paper will shed 
light on the nature of GRB progenitors. 

(4) In the Star Cluster Case, the EIC luminosity de¬ 
pends on the number density of massive stars at the time 
when the GRB occurs. There are large uncertainties in 
the first few M y r of the s tar clusters’ evolution (see e.g. 
IPortegies Zwart et al.ll 20 ld l. especially in the so-called Gas 
Expulsion Stage (GES) during which the residual gas in the 
cluster is blown away by stellar winds or outflows. The loss 
of mass during GES will cause the cluster to expand and the 
shallower potential well may allow a fraction of the stars (es¬ 
pecially the ones at outer radius) to escape. Therefore, the 
overall stellar number density drops with time. Long GRBs 
are concentrated on the very brightest regions of actively 
star-forming galaxie s, where young massi ve clusters are par¬ 
ticularly abundant jFruchter et all 12006!) . However, little is 
known about the very nearby (~ 1 pc) environment around 
the long GRB progenitors, and it is unclear whether the 
stars escaped 13 ! from the clusters. The EIC signal proposed 
in this paper could potentially give us some hint. For ex¬ 
ample, an upper limit of Leic < 10 49 erg/s will constrain 
the number density of massive (> 20 - 30A-/q) stars to be 
n» < 10 4 pc -3 . 

(5) There are a few potential issues in our numerical 
procedures developed in Section [2] (a) We assume electrons 
are accelerated to ultra-relativistic energies and have a sim¬ 
ple powerlaw Lorentz factor (LF) distribution. On one hand, 
if GRB prompt emission is produced via sub-photospheric 
Comtonization, then 7 m (or electrons’ thermal LF) might be 
of order a few or less. In that case, the EIC emission would 
be much weaker than the internal dissipation model con¬ 
sidered in this work. On the other hand, even if electrons 
are accelerated to a simple powerlaw, due to synchrotron, 
synchrotron-self-Compton, or EIC cooling, there’s a cool¬ 
ing LF 7 C above which electrons’ distribution turns softer. 
The resulting EIC spectrum at energy higher than r 7 c m e c 2 
will be softer than our prediction, but the total EIC energy 
will not change much, as long as 7 C > 7 m . (b) Electron’s 
LF distribution is actually evolving with time, i.e. 7 m , 7 C 
and possibly the powerlaw index p are all functions of time 
(or radius). I 11 principle, by convolving the time dependent 
LF distribution with single particle emission, we can calcu¬ 
late the EIC emission to a higher degree of accuracy. How- 


12 lHammer et all !2006t) found GRB 980425 (the nearest GRB 
at z = 0.008. e.g. iGalama et al.lll998l ) to be ~ 800 pc away from 
the nearby massive star-forming region. 
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ever, since we are not yet sure where the internal dissipation 
accelerates electrons an d how long will electrons stay hot 
ijKiimar fc Zhanell2015l h such a detailed calculation is left 
for possible future investigation. Instead, constant 7 m and 
p are used at all radius. As shown in Section EH (Fig GO), by 
using a constant 7™. at all radius, we are making an assump¬ 
tion that electrons are hot near the radius r e ff where most 
scatterings happen. 

(6) Another interesting application is the IC scattering 
of the Cosmic Microwave Background (CMB) at high red- 
sliift by Pop III GRB jets. The earliest Pop III stars form in 
minihalos of mass ~ 10® Mq at redsliift z ~ 20 - 30 (iBromml 
120138 1. The CMB photon number density is 

ncMB ~ aT 4 /(2.7kT) = 3.3 x 10 6 [(1 + z)/20] 3 cm“ 3 (50) 

where we have used the CMB temperature T = 2.7(1+ 2) = 
54 K (l + z)/20 (i.e. e 0 = 4.7x 10 -3 eV). The forward shock 
(FS) is highly relativistic at deceleration radius, and elec¬ 
trons are accelerated to LF 7 e = 1.8 x 10 5 (using t e = 0.1, 
and Tdec = 100). Therefore, EIC emission from the FS could 
be observed at 0.8 GeV F^ ec 2 . If we use a jet kinetic energy 
Ej = 10 55 erg (isotropic equivalent) and a uniform circum- 
burst medium density of n = 100 cm -3 , we get the decel¬ 
eration radius rd ec — 10 17 cm fEa 1371) . The total (isotropic) 
number of EIC photons from FS is 

n eic ^ {TFsncMB^nr 3 /3)\ r =r dea 

~ 3.0 x 10 52 [(1 + z)/20] 3 n 2 r 4 dec , 17 

Since N^f c strongly depends on r^c, a more accurate calcu¬ 
lation of the dynamical evolution is needed. Using N^fc = 
3 x 10 52 , we expect a photon fluence of ~ 5 x 10 -9 cm~ 2 at 
~ 1 GeV from redshift 20. 

(7) Lastly, we estimate the possible contribution of the 
EIC emission to the extragalactic 7-ray background (EGB) 
around 100 GeV, by the FS channel or ID channel when 
GRBs have large LF (e.g. T ^ 500). We optimistically es¬ 
timate the average (isotropic) EIC energy around 100 GeV 
from one GRB to be Eeic = 10 51 erg. The average (ob- 
ser ved) GR B r ate in the l ocal ( z < 1) Universe is estimated 
by IWanderman fc Piranl d201 (J l to be p ~ 3 Gpc~ 3 yr~ 1 . 
We integrate over the local Universe to Rmax = 4 Gpc 
and estimate the intensity of EGB contributed by the 
EIC channel to be Iegb — EEicpRmax/VK ~ 2 x 
10~ 6 MeVcm~ 2 s~ 1 sr^ 1 . The EGB at 100 GeV observed 
by Fermi LAT is ~ 10 -4 MeVcm~ 2 s~ 1 sr~ 1 I Abdo et al. l 
l2010l h Therefore, the EIC channel in GRBs could (at 
most) contribute a small fraction of EGB. The possible 
contr i butions from ot her GRB channels are estimated by 
ICasanova et al.1 d2007t l . 


5 CONCLUSIONS 

GRB progenitors may be surrounded by a significant exter¬ 
nal photon field (EPF), due to the existence of a massive 
companion star (Case I), a strong stellar wind (Case II), or 
a dense star cluster (Case III). We calculate the IC scatter¬ 
ing of the EPF (i.e. EIC emission) by the hot electrons in 
the GRB jet, which could be accelerated by internal dissipa¬ 
tion or external shocks. As long as electrons are accelerated 
to ultra-relativistic powerlaw, the results presented in this 


work is independent of the details of jet dissipation mecha¬ 
nism and electron acceleration process. 

In each Case, the EPF-contributing star(s) are assumed 
to have bolometric luminosity 10 39 erg/s and effective tem¬ 
perature 5 x 10 4 K. The GRB jet is assumed to have to¬ 
tal kinetic energy 10 54 erg (isotropic equivalent), duration 
Tj = 10 s and bulk Lorentz factor T. We assume that elec¬ 
trons are accelerated by internal dissipation to a powerlaw 
distribution with minimum Lorentz factor and index p. 
We present in this work the EIC lightcurves and spectra for 
each case, for a variety of EPF and jet parameters (F, 7 m 
and p). We take into account the equal-arrival time surface 
and possible absorption of high energy photons caused by 
pair production with the prompt 7-rays. 

For our fiducial jet parameters Y = 300, = 100 and 

p = 2.2, the EIC spectrum peaks at u p ~ 10 GeV. The EIC 
spectrum below v p is vL v oc v 2 and oc above v v \ pair 

production softens the spectrum at high energies (above the 
threshold for e 41 production). I11 Case (III), significant EIC 
emission also comes from electrons accelerated by external 
shocks, where we use e e = 0.1 as the fraction of internal 
energy shared by hot electrons in the shocked regions. 

In Case (I), using binary separations 10 14 —10 16 cm, we 
get EIC luminosities of 10 47 — 10 5 ° erg/s (the effect of pair 
production on high energy photon luminosity is included in 
this and all calculations below), peaking at separation 10 15 
cm. In Case (II), soft photons originally from the progen¬ 
itor star are scattered first by the electrons in the stellar 
wind and again by the jet. Using a W-R-star mass loss rate 
10 -5 Mq/pt and wind velocity 10 s cm/s, we get an EIC lu¬ 
minosity of 10 48 erg/s. I11 Case (III), we assume that stars 
in the cluster create a uniform EPF. Electrons in the jet 
are accelerated by both internal dissipation (ID) and exter¬ 
nal shocks. We include both the reverse shock and forward 
shock. Using stellar number density 10 4 pc~ 3 and cluster ra¬ 
dius 1 pc, we get EIC luminosities of 10 49 erg/s from ID, and 
10 49 erg/s from external reverse shock and 10 46 erg/s from 
external forward shock. The EIC emission from ID lasts for 
t 0 bs — Tj ~ 10 s, but that from external shocks lasts longer, 
approximately 200 s from reverse shock and 2 x 10 4 s from 
forward shock. 

We note that EIC emission from ID relies on the as¬ 
sumption that electrons are hot at the right radius r e ff. 
~ d (the binary separation, Case I), ~ rt r (~ 3 x 10 14 cm, 
the radius where the jet becomes transparent, Case II) and 
10 16 cm I'max, 16 (the maximum radius electrons stay hot, 
Case III). In addition, EIC emission from external shocks 
strongly depends on the dynamical evolution of the jet. 

Generally, from a GRB at redshift z = 1, a photon 
fluence of ICE 7 — 10 -4 cm~ 2 at 1 — 100 GeV is expected. 
Future observations (by e.g. Fermi LAT, HAWC and CTA) 
can put constraints on: (1) the existence of such EPFs as 
described in this paper and hence on the nature of GRB 
progenitors; (2) the radius where the jet internal dissipation 
process accelerates electrons. 
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APPENDIX A: 

There are two potential pair-production channels that could 
absorb the high energy (10 — 100 GeV) the EIC photons. 
(1) In the lab frame, the EPF (~ 10 eV) will interact 
with ~ 100 GeV photons. We show below that this chan¬ 
nel is usually not important. (2) In the jet comoving frame, 
the GRB prompt emission (~ 1 keV) will interact with 
~ 1 GeV photons (~ 100 GeV in t he lab frame). We 
consid e r this channel in detail, fo l lowin g iMadau fc Phimievl 
(l99dh : [Gilmore & Ramirez-Ruij li201Q j. and this channel is 
included in our numerical calculation (Section [2]). 
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A1 Absorption by EPF Itself 


The only concern here is for photons of energy ^ 100 GeV. 
The optical depth of pair production is 

T-y-y ~ n-yCrZ (Al) 


where l is the path length of the high energy 7 -ray through 
the EPF, n 7 is the number density of field photons near 
the energy of maximum cross section, and a ~ O.Ictt is the 
approximate cross section. 

I 11 the Binary Case, most scatterings happen at radius 
r ~ d, so we have 


(!) ~ Lb H %/2. 2 Lb ,39 

r 77 — I , jo \ -(l--)crcZ = 1.6x10 -- 

77 47r(d 2 + d 2 )eoc 2 di4£o,i 


(A2) 


In the Wind Case , most scatterings happen at radius r ~ rtr 
(see Eq |21l and !26l) . so we have 


r ( 2 ) 


LbT w i n d{rtr) 

47rr| r eoc 


artr = 2 x 10 6 


E2.5Eb,39A/-5 
Ej^ArjKN tO, 1 VS 


(A3) 


In the 5tar Cluster Case, EIC photons need to penetrate the 
whole cluster that has radius R = lpc R pc , so we have 


(3) n,LbR 


aR = 4.5 x 10 


eo c 


-i n*.4 Lf>, 39 
eo.i 


(A4) 


These are upper limits, since some EPF photons are either 
offset from the maximum cross section energy or at energy 
lower than the threshold. Anyway, the absorption is only 
important in Binary Case when d 10 12 cm and in Star 
Cluster Case when R 1 pc. Under these conditions, there 
will be a spectral cut off at 100 GeV. The absorption is 
not imp ortant in Wind Case. A mo re det ailed calculation is 
given bv lGilmore fe Ramirez-RuI5 (j201(t ). 


where q = -y/1 — u' th /v[. Assuming that low energy photons 
are moving isotropically in the jet comoving frame, we get 
the mean free path 

1 t 1 

,r) — - J d(cos6'i)(l - cosO'i) J , dvln^a^ (A 8 ) 

Then the pair-production optical depth is the comoving dy¬ 
namical timescale divided by the mean time interval between 
two collisions, i.e. 


TrrOV) 


r 

n 77 (z/', r) 


(A9) 


Therefore, only a fraction e -T77 of the EIC photons in 
Eg. (11311 can escape. A simple estimate of the pair-production 
optical depth is 


'T'y'y Csi 0.1 crr 


L-y /4 


4nr 2 chi/' th r 2 T 


- ~ 0 65 L ^,52{hv/10GeV) 


nsPl, 


(A10) 

where we assume 1/4 of all the prompt 7 -rays contribute to 
the absorption of high energy photons. We can see: (1) t 77 
is very sensitive to the bulk Lorentz factor (we have used 
T = 300T2.5); ( 2 ) most absorption happens at small radius 
(r 15 < 1) and high frequency (hu ^ 10 GeV). 


A2 Absorption by GRB Prompt Emission 


In the jet comoving frame, a high energy EIC photon of fre¬ 
quency v’ is surrounded by low energy photons (frequency 
v\), which ca n be described by the so-called Band function 
(iBand et al.lll993h . We simplify the Band function as a bro¬ 
ken powerlaw of indexes a\ and ct 2 , below and above the 
break frequency Vb- From the prompt emission bolometric 
isotropic luminosity L 7 , and the radius r where the low 
energy photons are produced, we get the comoving-frame 
specific number density 


n ',( r ) = _ Si _ ( _i_i_j" 1 

v i 4nr 2 chv£ an + 2 02 + 2 

r MM) ai if v'l < 

IMM)“ 2 if v'i > v' b 


(A5) 


where v' b = Vb/T is the break frequency in the comoving 
frame. The pair-production threshold is 


hv'th = 


2 (m e c 2 ) 2 
hv'(l — cos6[) 


(A 6 ) 


where 9[ is the angle between high(iz') and low(i/ i ') energy 
photons’ m omentum vectors. The cr oss section for pair pro¬ 
duction is (|Madau fc Phinnevlll996tf 


/ / / a /\ 3 (JT , , 2\ 

o--ry{v ,Vi,9i) = -(!-? ) 


(3 — qjln 


1 + g 
1-q 


+ 2 q(q 2 ~ 2) 

(AT) 



























